A small universe after all? 
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The cosmic microwave background radiation allows us to measure both the geometry and topology 
of the universe. It has been argued that the COBE-DMR data already rule out models that are 
multiply connected on scales smaller than the particle horizon. Here we show the opposite is true: 
compact (small) hyperbolic universes are favoured over their infinite counterparts. For a density 
parameter of Q = 0.3, the compact models are a better fit to COBE-DMR (relative likelihood 
~ 20) and the large-scale structure data (as increases by ~ 25%). 
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Measurements of the cosmic microwave background ra- 
diation (CMBR) provide a powerful probe of the geom- 
etry and topology of the universe. The geometry of the 
universe is reflected in the height, position and spacing of 
acoustic peaks in the CMBR angular power spectrum |Q , 
while the topology of the universe is betrayed by matched 
circles in the microwave sky g,|3|]. The topology of the 
universe also influences the angular power spectrum by 
"quantising" the spectrum of density fluctuations, and in 
many cases, by imposing a long wavelength cut-off. An 
infrared cut-off in the spectrum of density fluctuations 
translates into a suppression of large angle CMBR fluc- 
tuations on the surface of last scatter. This effect has 
been used to rule out flat [Q and hyperbolic || models 
with toroidal topologies. 

Earlier we conjectured |J that similar negative con- 
clusions would not apply to generic compact hyperbolic 
models as the majority of the large angle CMBR power 
in a hyperbolic universe comes not from the surface of 
last scatter, but from the decay of curvature perturba- 
tions along the line of sight. To a lesser extent, the same 
will be true in the newly popular flat models with a cos- 
mological constant, or in flat models with other forms of 
exotic dark matter such as tangled string networks [ JL8| . 
In what follows we shall not only verify our conjecture, 
but show that finite hyperbolic models actually provide 
a significantly better fit to the COBE-DMR data than 
their infinite counterparts. 

Our main focus will be on universes with compact hy- 
perbolic spatial sections || as these are the most appeal- 
ing from a theoretical standpoint. However, the growing 
body of observational evidence favouring a flat universe 
with a cosmological constant m prompts us to reconsider 
models with three-torus topology. We refer to a model 
as being "small" if its comoving spatial volume is less 
than the comoving volume enclosed by the particle hori- 
zon (measured in the covering space). The ratio of the 
horizon volume to the volume of the space exceeds 500 for 
several of the models we looked at. Small universe models 
are obtained from the usual FRW models by making iden- 
tifications between different points in space. These iden- 
tifications break global isotropy and homogeneity, but 
do not alter the evolution history. There is one caveat to 
the last statement: by altering the mode spectrum, the 
topological identifications will alter the vacuum struc- 
ture, leading to a Casimir-like vacuum energy that could 



alter the dynamics. If an effective cosmological constant 
could be linked to the universe having non-trivial topol- 
ogy, we would have a strong motivation for re-considering 
flat models. 

On large angular scales, temperature fluctuations in 
the CMBR are related to the fluctuations in the gauge- 
invariant gravitational potential <£ by the Sachs- Wolfe 
equation 
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Here 77 denotes the conformal time, $' = c^, and r] s i s 
and f] denote recombination and the present day respec- 
tively. The evolution of the gauge- invariant potential 
from last scatter until today is described by 

$" + 3W(1 + <? a )& - c 2 s V§ + (2W + (1 + c 2 s )H)§ = , 

(2) 

where H is the conformal Hubble factor and c s is the 
sound speed in the cosmological fluid. In a flat matter 
dominated universe we have c s = 0, H = 2/r] and (0) tells 
us that $' = 0. Consequently, the second term in ffl) van- 
ishes and the temperature fluctuations are all imprinted 
when matter and radiation decouple. The presence of ei- 
ther curvature or a cosmological constant alters the time 
evolution of the expansion rate H, leading to a decay of 
the potential <£. In these models, the line-of-sight integral 
(Integrated Sachs- Wolfe, or ISW effect) in ([I]) can be the 
dominant source of large angle temperature fluctuations. 
The gravitational potential can be expanded in terms 
of the eigenmodes, \£ g , of the Laplacian: 
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In the above equation n denotes the multiplicity of each 
eigenmode and 5 denotes the amplitude. The eigenmodes 
are found as solutions of the equation V\& q = —q 2 ^/ q . 
The geometry of the space determines the form of the 
Laplacian operator V, while the topology determines 
the boundary conditions. Though it is a simple exer- 
cise to write down the eigenmodes for any of the 10 flat 
topologies, hyperbolic manifolds have defied description. 



The first breakthrough came last year when Inoue found 
the 14 lowest eigenmodes of the Thurston space |l(| us- 
ing a numerical method (TTJ originally developed for 2- 
dimensional manifolds. By refining the method and using 
a more powerful computer, Inoue extended the count to 
include the first 36 eigenmodes. These modes have since 
been used to study the CMBR [[l2| in a universe with 
Thurston topology. Using the same numerical method, 
Aurich |Hj has studied the CMBR in a small hyperbolic 
orbifold with tetrahedral topology. Recently a new, fully 
automated, algorithm for finding the eigenmodes 111] 
was discovered and implemented. The list of solutions 
has grown from 36 to several thousand in the past two 
weeks. Our cosmological simulations use the first 100+ 
modes for each of 4 spaces selected from the SnapPea |lq| 
census of compact hyperbolic manifolds. Our selections 
are: the Weeks space, m003(-3,l), as it is smallest know; 
the Thurston space, m003(-2,3), in order to compare our 
results to Inoue's; and two larger examples, s718(l,l) 
and v3509(4,3), to see how the size of the space influ- 
ences the CMBR. In units of the curvature radius cubed, 
the volumes of our selections are 0.9427, 0.9814, 2.2726 
and 6.2392 respectively. To get a feel for how the eigen- 
modes look, Figure 1 shows the lowest eigenmode of the 
Weeks space extended across the entire Poincare ball. 
The Poincare representation is a conformal mapping of 
H 3 into a unit ball in E 3 . 




FIG. 1. The first eigenmode of the Weeks space. The three 
views in each panel are, respectively, the x = 0, y = and 
z — slices through the Poincare ball. The upper panel shows 
the corresponding cross sections through the fundamental cell. 

The cosmological simulations are readily performed in 
spherical coordinates, where the eigenmodes can be ex- 
panded: 



^q = Yl J2 A qirnRqt(r)Y em (6,<l>). 
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Here the Y^m's are spherical harmonics and the radial 
function R q £ are either spherical (flat space) or hyper- 
spherical (hyperbolic space) Bessel functions. It is worth 
emphasising how dramatically different the eigenmodes 



look in compact flat spaces and compact hyperbolic 
spaces. In flat space the expansion coefficients A q £ m 
are peaked around certain values of £, and we can al- 
ways rotate the coordinate system so that A q £ m = for 
m y^ 0. In other words, global isotropy is badly bro- 
ken. In contrast, the A q i m 's for hyperbolic eigenmodes 
are statistically independent of £ and m. The A q £ m 's are 
gaussian distributed pseudo-random numbers with vari- 
ance proportional to 1/k 2 |l4| . W hile it is possible to 
assign a wavenumber k = \/ q 2 — 1 to hyperbolic eigen- 
modes, it is impossible to assign a wavevector. They are 
essentially omni-directional. You do not need inflation 
to explain why the microwave sky is nearly isotropic in a 
small hyperbolic universe 0]. 

Once the eigenmodes have been found, it is a sim- 
ple though laborious task to evaluate equation (|l|) using 
equations ([|) and (||). We took the mode amplitudes to 
be gaussian random numbers with variance \fkj(k 2 + 4) 
(hyperbolic space) or 1/k 3 / 2 (flat space). These are the 
standard inflationary power spectra. 
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FIG. 2. The angular power spectrum AT = y/£(£+l)Ce/6 
for the Weeks universe (upper panel) and the corresponding 
infinite universe (lower panel). Both models have Q — 0.3. 
The shaded areas mark the 1-a cosmic variance intervals and 
the data points are from the 4-year COBE-DMR data. The 
dashed line in the upper panel is the surface of last scatter 
contribution to the angular power spectrum. 

Figure 2 shows the angular power spectrum for a uni- 
verse with density parameter Q = 0.3 and Weeks topol- 
ogy. Also shown is the angular power spectrum for 
the corresponding infinite model, along with the COBE- 
DMR data points |17J . The effect of varying the density 
parameter is shown in Figure 3, while the effect of vary- 
ing the volume of the space is shown in Figure 4. Notice 
that the large volume space v3509(4,3) produces an an- 
gular power spectrum very similar to the infinite model 
shown in Figure 1. 

We quantified how well each model fit the COBE-DMR 



data by modeling the spread in the AT's by smooth prob- 
ability distributions (the statistics are not gaussian) and 
performing a standard likelihood analysis. Our results 
are based on 1000 realizations of each topology for five 
values of the density parameter. The likelihoods relative 
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FIG. 3. Angular power spectra for the Thurston universe 
with Q = 0.2, 0.3, 0.4 and 0.5. These results compare well 
with those found in Ref. [12] 
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FIG. 4. Angular power spectra for 4 compact hyperbolic 

universes with Q Q — 0.3. 



TABLE I. Relative likelihoods. 
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to the corresponding infinite model are listed in Table 1, 
along with the likelihoods of the infinite hyperbolic mod- 
els relative to a fiducial flat matter dominated universe. 
The lower the density, the better the compact models 
fare relative to their infinite counterparts. The mecha- 
nism behind this result can be seen at work in Figure 1. 
As we decrease the density, the infrared cut-off in the 
mode spectrum reduces the contribution from the first 
term in equation (y) , while the line of sight contribution 
from the second term becomes increasingly important. 
Rather miraculously, the two effects almost precisely can- 
cel out for the low volume models, resulting in a flat or 
mildly tilted spectrum. On small scales there is no differ- 
ence in the shape of the angular power spectra for com- 
pact and infinite models, but the compact models have a 
higher COBE-DMR normalization. This helps raise the 
predicted size of density fluctuations on 8 h~ x Mpc scales 
from the low value of as = 0.6 for an infinite model, to 
the larger value of as = 0.75 for the Weeks model (both 
with Q = 0.3). Measurements of the present day cluster 
abundance (ll| favour as = 0.9 ± 0.1 if fl = 0.3. At this 
stage we should stress that our results primarily affect 
the fit to the COBE-DMR data. On small scales there 
is a growing body of observational evidence [p| for an 
acoustic peak at £ ~ 220, which is consistent with the 
universe being flat, not hyperbolic. We should also men- 
tion that our results disagree with those of Bond et al |2(]] 
based on the method of images, and agree with those of 
Inoue et al HI based on the finite element method. 
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FIG. 5. The microwave sky in a flat three torus universe, 
with (upper panel) and without (lower panel) a cosmological 
constant term. 



If the universe is flat, but dominated by a cosmological 
constant, could the toroidal universe models be saved by 
contributions from the the line of sight integral in (uh? 
Our preliminary investigation suggests that they may, 
but not so much in the power spectrum department as in 
the breaking of global isotropy. For example, a flat cubi- 



cal three-torus with side length equal to half the horizon 
distance has a relative likelihood of 0.26 when Q a = and 
0.25 when Qa = 0.7. The real difference between these 
two models can be seen in Figure 5. Because the ISW 
effect mixes together different modes sampled at different 
points, it helps to hide the egregious breaking of global 
isotropy that lead to the matter dominated versions of 
these models being ruled out Q. 

Returning to the compact hyperbolic models, we want 
to see if the ISW effect ruins the the matched circle 
test B for non-trivial topology. The matched "circles in 
the sky" occur wherever the surface of last scatter self- 
intersects. Since the surface of last scatter is a 2-sphere, 
the intersections occur along circles. We see two copies 
of each circle of intersection, centered at different points 
on the sky. The portion of the microwave temperature 
coming from the surface of last scatter will be identical 
around each circle. However, the ISW contribution will 
be uncorrelated. 



ever, the ISW effect only operates on large angular scales, 
so we filter Figure 6 to remove all power below £ = 21. 
The temperature match for the filtered sky is shown in 
Figure 7b. The correlation coefficient |J increases from 
0.29 to 0.92 after filtering out modes with £ < 20. The 
matched circle pairs will persist until the visible universe 
is simply connected, which occurs at around Q = 0.95 
for most models in the SnapPea census. If we do live in 
a small universe, the Microwave Anisotropy Probe will 
find matched circles in the sky when it starts collecting 
data in 2001. 

This work grew out of earlier collaborations and ex- 
tensive discussion with Glenn Starkman and Jeff Weeks. 
Financial support was provided by NASA through their 
fund ing of the Microwave Anisotrop y Probe satellite mis- 
sion [ittp : //map . gsf c . nasa . gov/ . 
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FIG. 6. The microwave sky in one realization of the Weeks 
universe with Q — 0.3. The white lines mark one pair of 
matched circles. 



a) o 




b) 



FIG. 7. The temperatures around a pair of matched circles 
before and after the low multipoles are removed. 

Taking one realization of the Weeks universe (Fig- 
ure 6), we find the temperatures around a pair of matched 
circles (Figure 7a) and see that the match is poor. How- 
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